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Thermodynamics of Lattice QCD with massless quarks and chiral 
4-fermion interactions.* 
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We have simulated lattice QCD with an irrelevant 4-fermion interaction and 2 zero mass quarks. The chiral 
phase transition is observed to be second order and we discuss extraction of critical exponents. 
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1. The action 

We have extended the standard staggered 
quark QCD action by the addition of a 4-fcrmion 
term with a £7(1) x ?7(1) chiral flavour symmetry, 
generated by (£5,75^5) @]. (For earlier, related 
work see [Q.) Introducing the auxiliary fields 
(cr, 7r) to render the action quadratic we have 

L = -/3^[1- ~Re(Tr D UT/W)] + -ly^ft A* Aij) 
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for our molecular dynamics lagrangian, where 
A=p + m+—^2(ai + ieTVi) (2) 
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The advantage of this action is that the Dirac 
operator remains non-singular when m — so 
that we can work in the chiral limit. Note that 
we have an exact chiral symmetry generated by 
75 £5 which means we will have a massless pion 
at to = when this symmetry is spontaneously 
broken. 

The 4-fermion term is an irrelevant operator, so 
this theory should have the same continuum limit 
— continuum QCD — as the standard action. 

2. Application to Nf = 2 thermodynamics 

Nf — 2 QCD is believed to have a second or- 
der finite temperature chiral phase transition for 
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m = which weakens to a rapid crossover, with 
no actual phase transition, for m 7^ 0. The crit- 
ical exponents that describe the neighbourhood 
of this critical point are expected to be those of 
an 0(4) spin model in 3 dimensions. (See |3| for 
a recent review of lattice QCD thermodynamics 
and critical exponents.) 

For the staggered lattice formulation the crit- 
ical exponents are expected to be 0(2) or 0(4). 
(For recent work on staggered critical exponents 
see jij.) Here the advantage of being able to work 
at zero quark mass is clear. We can hope to see 
clear evidence for the second order nature of the 
phase transition. In addition one can measure 
certain critical exponents directly. In particular 
one can measure the exponent (3 m from 



= const(P c - (3f™ = 1 (a) 



(3) 



where fi m = 0.5 for mean field theory, (3 m « 0.38 
for 0(4), and f3 m « 0.35 for 0(2). In the high 
temperature phase, the common n and a screen- 
ing masses increase from zero as 



const((3 — p c 
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where 7™ = 1 for mean field theory, 7 m « 1.44 
for 0(4) and j m w 1.32 for 0(2). The same criti- 
cal exponent governs the vanishing of the a mass 
below the transition. If we now fix (3 = (3 C and 
decrease the mass m from a finite value to zero, 
(ipip) vanishes as 

('ipip) = const m 1 / 5 (5) 

where 5 — 3 for mean field theory, S sa 4.8 for 
0(4), and 6 » 4.8 for 0(2). 
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Of the fluctuation quantities (susceptibilities), 
the most accessible is the "magnetic" susceptibil- 
ity 

x* = ^[<(£^))V<(£^))> 2 ] ( fi ) 

X X 

which scales as 

Xa = const(P - (3 C )~^ . (7) 

3. Lattice simulations and preliminary re- 
sults 

Our simulations are being carried out using 
hybrid molecular dynamics with noisy fermions 
to accommodate Nf = 2. Because of the 
exact flavour U(l) ax iai symmetry at to = 
0, the direction of symmetry breaking in the 
((iJip),i(^5^5fp)) space (or (ct,7t) space is not 
determined. On a finite lattice, this direc- 
tion rotates during the run, forcing us to use 
\J (tpi^) 2 — (V'TsCsV') 2 or \J ( (T ) 2 + (k) 2 as our or- 
der parameter on each configuration. (Here () 
should be taken to mean a lattice average, not 
an ensemble average.) This estimate differs from 
the true value by 0(1/ y/V), which can only be re- 
moved by working at more than 1 spatial volume 
for each N t . 

We are currently running on 8 3 x 4, 12 2 x 24 x 4, 
12 3 x 6, and 18 3 x 6 lattices at 7 = 20, and on a 
12 3 x 6 at 7 = 10. 

Whereas our earlier simulations at Nt = 4 
and 7 = 10 showed a first order transition, all 
Nt — 6 combinations above show evidence for the 
expected second order transition. (^V) an d the 
Wilson line show a sharp, but not discontinuous, 
transition with no sign of metastability, unlike the 
previous case. Figure |l| shows the (3 dependence 
of the chiral condensate for Nt — 6 and 7 = 20. In 
addition, close to the transition, observables show 
clear signs of critical slowing down with large fluc- 
tuations over many thousands of time units. (We 
have simulated for as many as 39,000 time units 
at a single (3 value.) Such a time evolution for the 
chiral condensate close to the transition is shown 
in figure ||. The Nt — 4, 7 = 20 runs show signs 
of critical fluctuations, but the transition occurs 
over a very narrow range and it will require more 



7 = 20 

0.5 1 




0.0 

5000 10000 15000 20000 

TIME 

Figure 2. Time evolution of ('tpip) at j3 = 5.42 on 
an 18 3 x 6 lattice at 7 = 20. 

statistics to tell if it has become second order or 
remains first order. 



3 



It is clear that we will need to extend our runs 
close to the phase transition. This is a conse- 
quence of the critical slowing down close to the 
second order transition. This currently limits our 
ability to extract accurate critical exponents due 
in part to the fact that our errors are poorly de- 
termined and partly because it prevents accurate 
extrapolation to infinite volume. In addition, the 
scaling window could be very narrow requiring us 
to consider more (i values close to the transition. 
Despite this the (V>V) data looks promising, and 
suggests a critical exponent (3 m ~ 0.3, however, 
this is on the basis of statistically poor fits. The 
N t = 6, 7 = 20 data shows relatively modest fi- 
nite volume effects. 

For N t = 4, 7 = 20, we find f3 c ~ 5.288. For 
N t = 6, 7 = 10, I3 C w 5.466, while for N t = 6, 
7 = 20, f3 c w 5.424. 

Figure |3| shows the a and 7r screening masses for 
the 18 3 x 6, 7 = 20 simulations. The a screen- 
ing mass shows the correct behaviour, decreas- 
ing to zero as (3 is increased through the transi- 
tion, and increasing from zero, degenerate with 
the 7r above the transition. Below the transition, 
the pion mass is consistent with zero. It appears 
doubtful that these masses can be determined ac- 
curately enough to obtain a good estimate for j m . 
The S mass can be calculated, but even assuming 
it remains distinct from the a/n mass above the 
transition, it will be difficult to tell how much 
of this is due to the explicit U{l) ax iai symmetry 
breaking provided by the 4-fermion interaction. 

4. Summary 

Nf = 2 lattice QCD with massless quarks 
and a weak 4-fermion interaction appears to have 
the expected second order transition, at least for 
N t > 6. More work is needed to clarify the N t = 4 
case. 

With more statistics the N t — 6 simulations 
should produce an accurate determination of the 
critical exponent f3 m . Moving to finite mass at 
(3 = (3 C should allow an accurate determination of 
S. 

Hadronic screening masses need further analy- 
sis. Other order parameters remain to be anal- 
ysed. 
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Figure 3. a and ir screening masses on an 18 3 x 6 
lattice at 7 = 20. 

Unfortunately, there is no obvious way to in- 
clude 4-fermion interactions with full SU(2) x 
SU (2) chiral flavour symmetry. 
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